Submitted 

Envelope solitons induced by high-order effects 
of hght-plasma interaction 

P. K. Shukla, R. FedeleEil, M. OnoratoS, N. L. TsintsadzeEl 

Institut fiir Theoretische Physik IV, FakuHdt fiir Physik und Astronomie, 
Ruhr-Universitdt Bochum, D-44780 Bochum, Germany 

The nonlinear coupling between the light beams and non-resonant ion density perturbations 
in a plasma is considered, taking into account the relativistic particle mass increase and the light 
beam ponderomotive force. A pair of equations comprising a nonlinear Schrodinger equation for 
the light beams and a driven (by the light beam pressure) ion-acoustic wave response is derived. 
It is shown that the stationary solutions of the nonlinear equations can be represented in the form 
^ ■ of a bright and dark/gray soliton for one-dimensional problem. We have also present a numerical 

' analysis which shows that our bright soliton solutions are stable exclusively for the values of the 

, parameters compatible with of our theory. 
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I. INTRODUCTION 



Recently, the investigations concerning the nonhnear dynamics governed by a multidimensional cubic-quintic non- 
linear Schrodinger equation (NLSE) received a great deal of attention. In this context, both localized vortex solitons 
. ^ ' ^^'^ nonlocalized optical vortices were studied Q. Cubic-quintic (2 + l)-dimensional NLSE were used to study the 
j • stability of spinning ring solitons and theoretical investigations to find solitary solutions for the cubic-quintic 
I " I (1 -|- l)-dimensional NLSE have been carried out. Dark solitary waves in the limit of small amplitude have been found, 
where the NLSE was reduced to a Korteweg-de Vries equation (KdVE) Moreover, both algebraic solitary wave 
^-H I solutions Q and traveling-wave solutions |^ have been found and criteria for existence and stability of soliton solu- 
K*" , tions have been established |^ . Additionally, a theory which connects envelope solitons of a wide class of generalized 
NLSEs with solitons of a wide class of generalized KdVE have been recently carried out for arbitrary amplitudes 
0; in particular, the theory was applied to find analytical bright, gray and dark envelope soliton solutions of the 
cubic-quintic NLSE and some other types of nonlinearities 0-0. 

It is well known that nonlinear interactions between intense laser beams and a plasma are responsible for numerous 
fS| nonlinear phenomena including parametric instabilities [ pO| , density cavitation, self- focusing and filamentation of light 
pl|-p^ , as well as the generation of wake fields iQ . Intense laser beams can cause density modifications through 
the ponderomotive force, enhance the electron mass due to relativistic effects and produce electron Joule heating. 
' ^ The interplay between the ponderomotive, relativistic and Joule heating nonlinearitries has been examined [Tl] ] in the 
(~| , context of laser plasma experiments and also in ionospheric modifications of the Earth's ionosphere by powerful radar 
beams. 

In this Letter, we investigate nonlinear interactions between the light beams and the non-resonant density pertur- 
bations, taking into account the combined effects of the light pressure induced ion density fluctuations and increased 
?-j [ electron mass. We show that, under suitable physical conditions for which our system can be described by a (1 -f 1)- 
dimensional cubic-quintic NLSE for the complex electromagnetic field amplitude, bright, gray and dark envelope 
solitons are analytically found. Finally, a stability analysis is carried out, which shows that our bright soliton solu- 
tions are stable. 



II. BASIC EQUATIONS 

We consider the propagation of a large amplitude circularly polarized electromagnetic wave with an electric field 
E = i?(x 4- iy) exp(—iLut + ikR), where uj is the wave frequency and k is the wavevector. The light equation in the 
presence of electron density perturbations in a plasma is obtained from 

An 1 

VB = — J + -dtF,, (1) 
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with 
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J = _e(no + ni)ve, 



(2) 



B = VA, (3) 

E = --dtA, (4) 
c 

and 

dtPe = -eE, (5) 

where B is the wave magnetic field, A is the vector potential, no and ni are the unperturbed and perturbed electron 
number densities, Vg is the particle quiver velocity induced by the photons, Pe = rngVe is the momentum, rUe = 
mo/(l — Vg/c^)^/^ is the mass, mo is the rest electron mass, e is the magnitude of the electron charge, and c is the 
speed of light in vacuum. The perturbation of the number density ni is reinforced by the light ponderomotive force. 
For our purposes, we have 

A 

(6) 



moC7e 

in view of (4) and (5). Here we have denoted 7e = y'l + e^A^/mQC*. Combining (1), (2), (3) and (6) we obtain 

A 
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d^A-c''V^A + wl{l + N) — = 0, (7) 



where Up = (47moe^/mo)^/^ is the plasma frequency, N = ni/riQ, and where we have introduced the Coulomb gauge 
VA = 0. We note that the term proportional to ni comes from the beating of the light quiver velocity with the 
electron density perturbations of the plasma slow motion, while the nonlinear term in "fe arises due to the electron 
mass increase in the light wave fields. Supposing that A = A5(r, r) exp(ikR — iuit) + complex conjugate, where r 
and T represent slowly varying space and time coordinates, we obtain from (7) 

2iu}{dr + VgVr)A, + c^VlA, + n^A, - ^ ^— = 0, (8) 

7e 

where = kc^/w is the group velocity of the light wave, and |9rA| << uiA has been invoked in view of the WKB 

approximation. We have denoted Q = We now derive the equation for low-phase velocity (in comparison 

with the electron thermal speed) density perturbations that are driven by the light wave ponderomotive force. The 
governing equations are the inertialess electron momentum equation 

= eVr?!> - moc2Vr7e - TeVr/n (rie/no), (9) 

the ion continuity equation 

drni + VriUiUi) =0, (10) 

and 

e Ti 

drUi + (ujVr)uj = Vr(p -V Jn {iii/no), (11) 

mi rrii 

where is the ambipolar potential, Uj is the fluid velocity associated with the plasma slow motion, and Te(Tj) is the 
electron (ion) temperature. The second term on the right-hand side of (9) represents the light pressure. Equations (9) 

to (11) form a closed system when the quasi-neutrality = Ui is invoked. The light ponderomotive force acting on 
the ion fluid is insignificant. Equation (9) shows that the electrons are pushed away from the region of maximum light 
intensity, and reinforce a space charge field (— V<?^) and the associated density fluctuations. The light ponderomotive 
force is transmitted to ions through the space charge field/ambipolar potential. Adding (9) and (11) and letting 
ne,i = no + Jii, Uj = Uio + = Uji, we obtain for m << no and | (uji-Vr) Ujij << I^T-Ujij 

drVia = -—C^Vrlel ' ^VrTli, (12) 

m, no 
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where for consistency we have supposed c^A^/ttiqC'* << 1 and, consequently, introduced a small perturbation 7ei of 
the electron relativistic factor 7e (7e « 1 + 7ei, where 7ei = e^Ag/2TOQC^) and Cg = [{Tg + Ti)/mj\^^'^ is the effective 
sound speed. According to the above approximations, Equation (8) becomes 

2iu:{dr + VgVr)A, + c^VlAs + - Lol) As - i^l [{N - 7el) - iV7el] A, = 0, (13) 

and (12), combined with the linearized version of (10), yields 

dlN - ClVlN = ^c'Vhei. (14) 

Equations (13) and (14) are the desired equations for coherent light beams that are coupled with non-resonant density 
perturbations in an clcctron-ion plasma. Note that, in principle, the quantities {N — j^i) and A^7ei, involved in Eq. 
(13), could be of the same order. In the following, this physical circumstance will be considered and, to this end, we 
seek possible stationary nonlinear solutions of equations (13) and (14) . in the form of envelope solitons 



III. ENVELOPE SOLITONS 



We introduce £, = r — Vr, where V is the velocity of the nonlinear waves, and assume A^ = a(^) exp(— zfior), where 
Oo is a constant. Hence, we readily obtain from (13) and (14) 

2iLo [(-V + Vg) -V^] a + c'V^a + (fl^ _ ^2 ^ 2u;l]o) a - [{N - 7ei) - N-fei] a = 0, (15) 

and 

(V-Ve)' N - Cl V|7V = V|7ei. (16) 

Tfli 

For the sake of simplicity, we consider here the one dimensional case for which can write {Y-V^)^ N = V^d^N. 
Consequently, (16) can be immediately integrated, giving 

rrii [V^ - Cj) 

By taking into account the explicit expression of 7ei, choosing |V| = |Vg| >> Cg, and combining (15) and (17), we 
easily get 

^dl^ + - + (;2|*|4] * = 0, (18) 

where we have introduced the following dimensionless quantities 

r/ = c ^/ujp, 



* = 



e a 



and 



92 
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with /z = rriQC^ /mi{V'^ — C^) « rriQC^ /rmV'^. 

Note that, since > 0, 92 is a negative quantity. Additionally, since N and 7ei must be of the same order, from (17) 

is evident that is of the order of the unity ^ 1), but slightly greater than 1. This circumstance is satisfied when 

1/2 

we choose, consistently, a group velocity Vg = V ^ {mQ/nii) c. This justifies why we have kept both the nonlinear 

terms in (18); accordingly, the terms gil^P and g2|*|^ are of the same order. In particular, if fi is exactly equal to 

1/2 

1 (i.e., we have exactly V = {mo/rrii) ' c), the (18) becomes 



1 



<92* + A2*-92|*|** = 0, (19) 



which shows that a part of the relativistic mass variation nonlinearity is exactly cancelled out by the light pondero- 
motive force driven supersonic electron density contribution. 
If we put 

*(?7,s) = *(r/)exp(-iA^s) , (20) 
where s is a new dimensionless timelike variable, the (18) can be cast as 

za,* + ia^* - + g2|*|^] * = 0, (21) 

Let us suppose that fi is (slightly) larger than the unity. In this way qi > , (21) admits bright, gray and dark 
envelope solitonlike solutions. In fact, from the results of recent investigations 0] that have found a wide class of 
envelope solitonlike solutions of (21), one can find, through (20), the following solitonlike solution for \l/(7y) 



^'(77) = y^u [l + e sech (77/ A)] exp l^iB 



, . 2. ^^^^^y(.-l)tanh(,/2A) 



+ icj)o} , (22) 



where (/)q is an arbitrary real constant, u — — 3gi/(8g2) = 3(/i — l)/(8/i), the constants e, A and B are given, 
respectively, by 



3gf V 3{p-i) 
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A = 2a 2 



64 192 1 



1/2 



2J2 



3(a^-1)^ 
128 \n\ 



1/2 



B ^VoA , 

provided that 

2 ^ 15(m-1)^ Vq^ 

128|/i| 2 

and the real constant Vq satisfies the condition 



/3(,-1)2 



64|Ai| V 64|Ai| 



Accordingly the definition of A2 implies that 



_ c2fc2 + 2ojno = ( 1 + + Fo' 

\ 64/1 

which is a condition for the real constant Qq. According to terminology and the results of Ref.s 0J^, we can distin- 
guish the following four cases, (a). < e < 1 (Vq 7^ 0): up-shifted bright soliton. 
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u{i] = 0) = u(l + e) , and lim u{ri) — u 

rf — >±oo 

which corresponds to a bright sohton of maximum amphtude (l+e)u and up-shifted by the quantity u. (b). — 1 < e < 
(Vb 7^ 0): gray soliton. 

u{ri = 0) = u(l — e) , and hm u{rf) — u 

7] — >±00 

which is a dark sohton with minimum amphtude (1 — e)u and reaching asimptoticaUy the upper hmit u. (c). e — 1 
(Vb = 0): upper-shifted bright soliton. 

u{ri = 0) = 2u , and hm u(ri) = u 

which corresponds to a bright sohton of maximum amphtude 2u and u p-shifted by the maximum quantity u. (d). 
e — —1 (Vb — 0): standard dark soliton. 

u{ri = 0) = , and hm u{ri) = u 

ri — ^±oo 

which is a dark sohton (zero minimum amphtude), reaching asimptoticaUy the upper hmit u. 
Correspondingly, the (17) gives the foUowing sohtonhke solution for the density fluctuation 



iV(ry) = M l*(^)l' 



(Ai-l)[l + esech(77/A)] 



On the other hand, according to Ref.s 



the (19) has the following bright envelope sohtonhke solution: 



6\Eo\ 



1/4 



sech 



1/2 



\/2\Eo\r] exp(i0o) 



(23) 



(24) 



where (j>o is an arbitrary real constant and Eq is a negative real constant satisfiying the condition = Eq . The 
latter implies the following condition for Hq uJ^ — c^fc^ + 2a;f2o + (2|£'o| — 1) = . Furthermore, the (17) implies 
that now the density fluctuation corresponds to the following sohtonhke solution 



7V(r7) = p. \^(vi)\^ = 



6\Eo\ 



1/2 



sech 



V2\Eo\v 



(25) 



We now investigate the stability of plane wave solutions of the one-dimensional equation (21). We first note that if we 
set q2 = 0, the corresponding equation is the well known defocusing cubic Schroedinger equation which is known to 
be stable. It is therefore interesting to study if the 92! 'I' I** term can modify the instability. The analysis is performed 
by seeking a solution corresponding to a uniform wave train perturbed by small disturbances: 



$ = [$o + /o(s,7y)] expi{hgi|$o|' - q2\^o\%s + e{s, t])}, 



(26) 



where p and 6 are considered to be small amplitude and small phase perturbations. We then substitute the perturbed 
solution in (21) and retain only the linear terms in p and 9. Since the resulting equation is linear we can now assume 
a solution for the perturbation of the form p = po expi[Kri — fls] and 6 = 60 exj)i[Kj] — fis]. The resulting dispersion 
relation is the following: 



f^' = ^(i^'+4(Zi|$o|' + 8g2|$o|') 



(27) 



This shows that the wave train is unstable if the perturbation K lies in the range of < -ft' < 2|<I>o|\/^'?i ^ 2g2|'i'op- 
According to the definition of gi and q2, the instability will occur only if/i > 1/(1 — 2|$oP)- The maximum instability 
occurs at if = \^o\^/^2q^^^Iq^\^^. 
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IV. NUMERICAL SIMULATIONS 



In this section we analyse numerically both the influence of the quintic nonlinearity on the modulational instability 
and on the stability of a class of solitonhke solutions obtained in the previous sections. Equation (21) is solved 
numerically using a standard pseudo-spectral code with a second order Runge-Kutta method for advancing in time. 
We recall that the use of pseudo-spectral code implies the assumption of periodic boundary conditions. 



A. Modulational instability 

Accordingly to the linear stability analysis performed previously, initial conditions for our numerical simulations 
are given as follows: 

<i>{x,0) ^(!>o[l + sCos{Lx)], (28) 

where £ is the amplitude of the small perturbation and is taken as 10~^ the amplitude of the umperturbed wave. 
Without loss of generality in our simulations we have choosen <I>o = 1 and L = 1 and have considered only one period 
of the perturbation. We have performed several numerical simulations with different values of the parameter /i. For 
$0 = 1 the theory predicts stability for /i < — 1. In Fig. 1 we show the evolution of a plane wave in the 77 — s plane for 
/i — —1.5. The initial wave field persists for all times. For ^ > 1 modulational instability should occur. In Fig. 2 we 
show the case of /i = —0.2; analogously with the standard modulational instability we observe a Fermi-Pasta-Ulam 
recurrence: periodically the perturbation grows, the wave reaches a maximum amplitude and then goes back to the 
initial condition. For /i > a completely different physics takes place: in Fig. 3, obatined for /x = 0.5, we do not 
observe anymore a recurrence and as time passes the wave amplitude increases while the its width decreses. This 
phenomenon corresponds to the initial stage of a wave collapse (see llH]). 



B. Stability of solitonlike solutions 

It is well know that the cubic NLS equations (92 = 0) posses solitons solutions if qi is greater than zero. We 
here investigate numerically if the solitonlike solutions described previously are stable or not. In order to to that we 
simply consider a solitonlike solution at time s = and we let evolve numerically the eq. (21). For semplicity we 
restrict our analysis to a sub class of solutions which corresponds to the case of bright solitons with Vq — 0. We have 
performed many different numerical simulations with different values of the parameter fj.. The major result obtained 
is the following: if /i > 1 solutions are stable and for < ^ < 1 are unstable. This is due to the fact that for 
yU > 1 the coefficients in front of the cubic and quintic nonlinearities, respectivelly qi and (72, have opposite sign and 
therefore there is a sort of balance between nonlinearities that stabilizes the solitonlike solution. This is not the case 
for < < 1: both nonlinearities have the same sign and the dispersion is not strong enough to balance them. In 
Figs. 4 and 5 we give numerical evidence of the results presented for /i = 1.1 and for fi = 0.9 respectivelly. The first 
case, Fig. 4, corresponds to a stable solutions (the wave profile does not change as time s passes). The second case. 
Fig. 5, is the unstable case: a clear increase in the wave amplitude is shown. 

According to the above stability analysis we can conclude that our soliton solutions are stable in the range where 



/i > 1 only; this inequality, according to section III, is consistent with the conditions for their existence in our problem. 



V. CONCLUSIONS 



To summarize, we have considered the nonlinear interaction between intense light beams and the non-resonant 
density perturbations, taking into account the relativistic mass increase of the electrons as well as the light beam 
ponderomotive force that reinforces the density perturbations in an electron-ion plasma. The nonlinear coupling is 
governed by a pair of equations which, in one space dimension admit stationary solutions in the form of a planar 
bright and dark/gray envelope solitons. The condition of stability of the bright soliton-like solutions has been found 
numerically and it has shown that they are stable just for the range of parameters required in our problem. 
Unfortunately, the numerical stability analysis for the family of dark and gray solitons, although revealed to be much 
more difficult, is on the way. Authors expect that the preliminary numerical analysis of the stability carried out in 
the present paper will be extended to the dark and gray solitons in a future work. 
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